arXiv:1506.00942vl [astro-ph.GA] 2Jun2015 


Local stability of self-gravitating disks in f{R) gravity 

Mahmood Roshan* • Shahram Abbassi^’^ 


Abstract In the framework of metric f{R) gravity, we 
find the dispersion relation for the propagation of tightly 
wound spiral density waves in the surface of rotating, self- 
gravitating disks. Also, new Toomre-like stability criteria 
for differentially rotating disks has been derived for both 
fluid and stellar disks. 
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1 Introduction 

The stars, gas and dust clouds in galactic disks congregate 
in spiral patterns and make bountiful structures in the uni¬ 
verse. Despite spiral arms beauties and many decades of 
concerted scientific investigation, much about them remain 
mysterious. Although the precise theory explaining the ori¬ 
gin and evolution of the spiral structures in spiral galaxies 
is not fully understood, it is widely agreed in the relevant 
literature that these patterns are gravitationally driven den 


sity waves in t he stellar disks (IBinnev & Tremainel 12008 


Sellwoodll2014h . Therefore, the density wave theory is an 
important tool for studying the dynamics and the evolution 
of spiral galaxies. In the frame work of this theory and 
using some approximations and assumptions. Alar Toomre 
showed that the differentially rotating disks in Newtonian 
dynamics is stable to all local axis ymmetric distu rbances if 
the dimensionless quantity Q > 1 (iToomrell 19641) . where Q 


is the Toomre’s Q parameter (see equations (l22li and (|2^ 
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for the definition of this parameter for fluid and stellar disks 
respectively). When Q < 1, on the other hand, thermal 
pressure and rotation are unable to stop the collapse of over- 
dense regions. 

In fact, he assumed that the density waves are tightly 
wound and the stellar orbits of the unperturbed disk are 
nearly circular. Despite these restrictive assumptions used 
in deriving this result, it has proved to be remarkably true 
and widely applicable. However, numerical studies of disk 
galaxies reveals that this criterion is not enough for complete 
stability of the stellar self-gravitating di sks (lHohllll97lf 


Miller et ^ 197(]l : Ostriker & Peeblesll973 ). In other words, 
these N-body simulations showed that Toomre’s criterion 
does not provide the global stability of the stellar disks and 
only guarantees the local stability of them. More specifi¬ 
cally, it turned out that simple models of rotationally domi¬ 
nated stellar disks are globally unstable to a pressure domi¬ 
nated bar-like structure. 

This kind of gravitational instability can be directly 
linked to the dark matter problem in spiral galaxies. This 
link has been known since 1973 when Ostriker and Peebles 
pointed out that a dark matter halo surrounded the galaxy 
may be an important sta bilizing agent of galactic disks 
(lOstriker & Peebleslll973[) and preyent the rapid bar forma¬ 
tion. Howeyer, it should be noted that existence of a spher¬ 
ical halo around a rotationally dominated (or cold) disk is 
not the only way to stabilize the disk. There are ot her possi¬ 
bilitie s which can proyide the global stability, see (ISellwood 
20141) for a reyiew of the subject. For example, e xistence of 
a massiye central bulge would stabilize the disk (Sellwood 


1985llSellwood & Moorel[l9^ISellwood & EyandboOlh. 

The main aim of this paper is to And the general¬ 
ized Toomre’s local stability criterion in the context of 
metric /(i?) grayity. As an example of such a study 
in other modified grayity theories see ( MilgromI 19891) 
where the Toomre’s criterion has been deriyed in the con- 
text of modified Newto nian dynamics (MOND). Also, 
Roshan & Abbassil (12015h haye deriyed the Toomre’s stabil- 
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ity criterion in the context of MOG. It is worth to mention 
that MOG is a Scalar-Tensor-Vector theory of gravity pre - 
sented to address the dark matter problem ( Moffat 2006h . 
Although a considerable amount of work has gone into sta¬ 
bility criteria, no study has been performed for investigat¬ 
ing these criteria in the context of f{R) gravity. However, 
it should be noted that the dynamical stability of spherical 
systems, and also the stability of spherically symmetric so- 
lutions in f(R) gravity h ave been widely in vestigated (see 
Noureen & Zubaiii (120151) and lSeifera (l2007h and references 
therein for more details.) Certainly, the results of this pa¬ 
per could be useful in the numerical simulations of the disk 
galaxies in f{R) gravity. 

It is worth mentioning that metric f{R) gravity is one 
of the simplest modifications to Einstein’s General Relativ¬ 
ity (GR). Strictly speaking, the generic action of this theory 
can be simply obtained just by replacing the Ricci scalar R 
with an arbitrary general function of R, namely f{R), in the 
Einstein-Hilbert action. Several aspect of this theory have 
been investigated in the literature. Interests to this theory 
increased when Carroll et al propo sed that f(fl ) gravi ty can 
solve the cosmic speedup enigma Carroll et^ (12004 ). Al¬ 
though this theory is known as a dark energy model, it has 
been applied to address t he dark matter pro b lem in the galac¬ 
tic scales. Eor example ICapozziello et all (120071) applied a 
power-law f{R) gravity to explain the rotation curves of 
low surface b rightness galaxies. We refer the reader to re- 
view papers dSotiriou & Earaorii 2010l : Clifton et al 2012 ) 
for more detail about this theory and its status among the 
other extended theories of dark energy. 

The structure of this paper is as follows. The weak field 
approximation of metric f{R) gravity is reviewed briefly in 
section |2] Also, the coupled Boltzmann and modified Pois¬ 
son equations are derived. In section [2 by linearizing the 
held equations, we derive the dispersion relation of tightly 
wound spiral density waves in both fluid and stellar disks. 
Eurthermore, using the new dispersion relations, we derive 
the local stability criterion. In fact, we And the generalized 
version of the Toomre’s stability criterion in the context of 
metric f{R) gravity. This new Toomre’s criterion is the 
main result of this paper. Einally, in section 0] results are 
discussed. 


where R,G, g, Sm are the Ricci scalar, Newton’s gravita¬ 
tional constant, the determinant of the metric tensor and the 
matter action, respectively. Note that throughout this paper 
we work in the system of units where the speed of light is 
c = 1. Eurthermore, we will assume that f{R) is analytic 
about i? = 0 so that it can be expanded as a power series 






( 2 ) 


There are some f{R) models of cosmological interest in the 
literature that can be expressed in such a form, for example, 
Starobinsky’s model Starobinskv (2007) 


f{R) = R + PRo 


can be expressed as 


1 + I “ 1 

-fin 


f{R) = R--fR 

Uq 


n2 , Prn{m + 1) , 


2Rl 


-R‘^ 


(3) 


(4) 


Variation of action ([T]) with respect to metric yield to the 
metric field equation: 


f'R^, - = SttGT^, 


(5) 


where prime denotes derivative with respect to R, i.e. 
f'{R) = and is the energy-momentum tensor. 
Throughout this paper, we restrict ourselves to the adia¬ 
batic approximation in which the evolution of the universe 
is very slow in comparison with local dynamics. This as¬ 
sumption allows us to use the Minkowski space-time instead 
of more natural choices such as the Eriedmann-Robertson- 
Walker (ERW) spa ce-time as the background metric (see 
((Clifton et alll2012h and references cited therein). It should 
be noted that in the current paper we consider the local grav¬ 
itational stability of the self-gravitating disks. This means 
that only the local dynamics of these systems matters for us. 

Therefore, in order to write the held equations in the 
first perturbation order, we decompose metric into the 
Minkowski metric = diag(-|-l, —1, —1, —1) plus a 
small perturbation h^i, \ <C 1) as follows 


2 Weak field limit of f{R) gravity 

Here, we briefly review the weak filed limit of f{R) grav¬ 
ity theory. The weak field limit of this theory has been 
widely i nvestigated, for a comprehensive review of the sub 
ject see (ISotiriou & Earaonill2010tlNoiiri & Odintsovll201 1 


De Eelice & Tsuiikawal ioiot Capozziello & De Laurentis 

201 ih . The general action for this theory is given by 


S = 


IGttG 


J V^f{R)d'^x + Sm 


( 1 ) 


= (1 + ~ (1 ~ 2'^){dx‘^ + dy'^ + dz^) ( 6 ) 

where doo = 2$(x, t) and hij = 2T'(x, f)(5y . We use the 
transverse gauge, i.e. = 0. where the traceless tensor 
Sij is known as the strain and is given by 




(7) 


Einally, substituting (|6]l into field equation (0 and keeping 
only terms linear in perturbations $ and fl', one can find 
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the following equation for the (0,0) component of the field 
equation 

/''V4(2^' - $) + = SttGp (8) 

where /q"^ = ^^\r=o and Similarly, 

preserving only the first-order perturbations, the trace of the 
field equation (|5]l takes the following form 


2.1 The modified Poisson equation in /(i?) gravity 

In this section we derive the Poisson equation in the weak 
field limit of metric f{R) gravity. For this purpose, it is con¬ 
venient to combine equations (H) and (|9]l as a single equa¬ 
tion. The result is the analouge the Poisson equation in 
Newtonian gravity. Using equations (|8ll and (|9|l, one can 
show that 


3/oV'‘(2T'-$)-f/oV^(2T'-$) =87rGp (9) 


It must be mentioned that we have assumed the perfect 
fluid’s energy-momentum tensor for T^,j. Thus p is the rest 
frame matter density. Equations (|8]l and @ all together with 
the continuity equation and the Euler, make a complete set 
of equations governing the dynamics of a fluid system in 
the weak filed limit of f{R) gravity. In fact, equations ® 
and (|9]l are the analog of the Poisson equation in Newtonian 
gravity. It is obvious that we have to And the corresponding 
equations for continuity equation and the Euler equation in 
the frame work of f{R) gravity. To do so, it is needed to 
take into account that metric f{R) theory is a metric the¬ 
ory of gravity. It is known that in metric theories of grav¬ 
ity the conservation of the energy momentum-tensor holds, 
i.e. V = 0. Consequently, non-spinning test parti¬ 
cles move on the geodesics of the metric tensor There¬ 
fore, the equation of motion of a test particle in this theory 
is given by 


dx^ 


( 10 ) 


where r is the proper time along the world line of the test 
particle. Using the metric component (|6]l, it is straightfor¬ 
ward to write Eq. (flOl i to Arst order in perturbations 


d^r 


= -Vfl> 


( 11 ) 


thus, only the metric potential $ directly appears in the 
equation of motion of a test particle. In other words, only 
$ appears in the gravitational potential. This is the case 
also for the continuity and Euler equation. Writing the di¬ 
vergence of the energy-momentum tensor, = 0, in 

the Newtonian limit, one can easily verify that the continuity 
and Euler equations are 

dp 

^ + V • (pv) = 0 (12) 


^2.^ SttG 
= -V^T> -f —^p 

Jo 

Now, substituting equation (fT4l i into @, we And 


(V^ — mQ)V^$ = —47rGa ( mpp — -V^p 


(14) 


(15) 


where a = 1//q and tuq is defined as 


A 

3/^' 


(16) 


in fact, equation (fTSl l is the modified version of the Poisson 
equation. Therefore, loosely speaking, we are dealing with a 
theory with two free parameters a and /3. In the limit uiq 
oo and a I, the standard Poisson equation is recovered, 
i.e. equation (fTSl l takes the form = AnGp. Solution of 
this fourth order equation (fTSl) can be written as follows 


$(r) = -4^Ga / / Gi(r",r')G 2 (r,r") 


X ^mop(r') - d^v'd^v" 


(17) 


where Gi and Gq are the Green functions; 


Gi(r",r') = - 
Gi(r,r") = - 


1 

47r |r' — r"| 

1 1 


47r |r — r" 


(18) 


it is natural to expect that the free parameters a and mo 
should be Axed by relevant observations. Also, in order to 
prevent the oscillatory solutions, mp should be positive, i.e. 

mp > 0. 

It is worth noting that the gravitational potential of a point 
source takes the following form 


$(r) = 


GMa 

r 


(19) 


Ov 1 

—-f(v-V)v = —Vp-Vfl> (13) 

ot p 

Equations (I8]l,(l9|l, (fT2l) and ( fOl l together with the equation 
of state relating p and p, make a set of equations governing 
the dynamics of a self-gravitating fluid system. 


This result can be straightforwardly derived from equa¬ 
tion dnii just by substituting the matter density as p(r) = 
MS{r), where M is the mass of the point particle. 
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2.2 The collisionless Boltzmann equation in f{R) gravity 

In this section, we derive the collisionless Boltzamnn equa¬ 
tion in the weak field limit of metric f{R) gravity. This 
equation is the main equation governing the dynamics of a 
stellar system. Assume that f{x^, u^) is the particle’s phase 
space distribution function. Where is the four ve¬ 

locity of the particle. In Newtonian dynamics the Boltzmann 
equation can be expressed as L\f] = 0, where L is the Liou- 
ville operator. It is worth remembering that in metric f{R) 
gravity non-spinning test particles move on the geodesics. 
Therefore with the aid of the geodesic equation, we obtain 
the covariant, relativistic generalization of the Boltzmann 
equation 

Lm = = 0 ( 20 ) 

It is straightforward to linearize this equation using the per¬ 
turbation introduced in Eq. (ISl). Keeping in mind that in the 
first order Newtonian limit Tqq ~ ^ and u* ~ u* = 
we have 

df df 

^+ vV /- V $-^=0 ( 21 ) 

as expected, only metric potential $ appears in the Boltz¬ 
mann equation. Mathematically, this equation is the same 
as the Boltzmann equation in Newtonian dynamics. How¬ 
ever, we know that the field equation governing <I> is differ¬ 
ent from the standard Poisson equation. 


3 Local stability of differentially rotating disks in f{R) 
gravity 


As we mentioned in the introduction, our main purpose in 
this paper is to find a Toomre-like stability criterion in the 
context of metric f{R) gravity. We remind the reader that 
Toomre’s criterion for a fluid (gaseous) disk is expressed as 


Qg 


VsK 

ttGE 


> 1 


( 22 ) 


where Vs is the sound speed in the fluid, k is the epicycle 
frequency and E is the surface matter density of the disk. 
Similarly, in the case of a stellar disk this criterion can be 
written as 


Qs 


3.36GE 


> 1 


(23) 


in this case a is the radial velocity dispersion. In Newto¬ 
nian gravity, Toomre’s criterion has been derived for differ¬ 
entially rotating disk where the tight-winding or the WKB 
approximation is satisfied. In fact, the tight-winding (or 
WKB) approximation provides some simplifications in the 
story and, without much loss of generality of the problem. 


makes analytic description possible and much simpler. In 
this approximation, the density wave in the disk can locally 
be regarded as a plane wave. This simplification enables 
us to find the dispersion relation for this kind of density 
waves and to describe their propagation in the disk. There¬ 
fore, WKB approximation is an indispensable tool for un¬ 
derstanding the behavior of density waves in rotating disks 
in the context of Newtonian gravity. As we shall see in 
the next sections, in order to do complete stability analy¬ 
sis in f{R) gravity and find the generalized version of the 
Toomre’s criterion, we also need to benefit this approxima¬ 
tion. 


3.1 Self-gravitating fluid disk 


In this section, we analyze the behavior of a tightly wound 
density wave in the framework of metric f{R) gravity. First, 
we use the modified Poisson equation (fTSl l in order to cal¬ 
culate the gravitational potential of a tightly wound surface 
density. Furthermore, we linearize the relevant equations in 
order to find the dispersion relation for such a density wave. 

The continuity equation (fTSI) . the Euler equation ( fT3] l and 
the modified Poisson equation ( fTSl l are the required equation 
for studying the gravitational stability of the system. It is im¬ 
portant to remember that we also need the equation of state 
of the fluid to make a complete description. We assume that 
the background disk is axisymmetric and is placed in 2 = 0 
plane. Using the cylindrical coordinates {R, cj), z), the con¬ 
tinuity equation reads (note that hereafter i? is a coordinate 
and should not be confused with the Ricci scalar) 

By i B i B 

^ + + (24, 

where E is the surface density and vr and are the velocity 
components in the radial and azimuthal directions respec¬ 
tively. Since the disk is located in the x — y plane, there are 
only two components for the Euler equation. These compo¬ 
nents can be , respectively, expressed as: 


dvR Bvr Bvr vl 

dt BR ^ R 3^ R 

8v^ Bv^ Bv^ v^vr _ 
^^ 'r~^ ^ R ~ 


-A($. 

BR '' 


h) (25) 


iA 

R'^ 


($ -h /i) (26) 


where we have chosen a simple barotropic equation of state 
as p = KTB, where K and 7 are constant real parameters. 
Also, h is the specific enthalpy defined as: 



(27) 


Now, let us perform the following perturbations: E = Eq + 

El, Ufl; = Vrq -b VRI = VRI, Vtj, = V^o V^i, $ = T>o -b ^>1 

and h = hQ-\- hi. We denote the equilibrium quantities with 
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a ”0” subscript and the perturbations with a ”1” subscript. Regarding the form of Si and $i, any solution of equa- 

After linearization, the fluid equations (l24l i and (l26l l can be tions (I28l l- (l30l l can be expressed as 
written as : 


9Si 


1 f) r)'y 


So dv^i 

'r~^ 


VRi{R,cl>,t) = VRa{R)ed--^-^*^ 

(37) 

0 (28) (i?, t) = 

(38) 

4 (E,(/),f) = 4(E)e*('"^--‘) 

(39) 


dvRi 

dt 


+ n 


dVRl 

d(j) 


— 2VtV^l — (‘^*1 + ^l) 


dv^i 


dcj) 




(29) In Newtonian gravity, using the properties of the WKB ap¬ 
proximation, equations (l28]) - (l30l) take the following forms 
respectively (for more detail see Binney & Tremaine 2008 ) 

(mfl - u;)S„ + fcSoUfia = 0 (40) 


where Vl{R) is the circular frequency and the epicyclic fre¬ 
quency K is 


4(E) 




(31) 


Also, the modified Poisson equation can be linearized as 


(V^ - = - 47rGct 


X 




:V2S4(z) 


(32) 


Now, consider an arbitrary spiral surface density perturba¬ 
tion mode as follows 




(33) 


where uj is the frequency of the mode, E is a slowly varying 
function of R and F{R) is the shape function. Also, m is a 
positive integer, and the perturbation has m-fold rotational 
symmetry. The WKB approximation requires that ^ 
1, where k{R) = ^- In this approximation the pitch angle 
of the spiral arms is very small at any radius. Therefore, one 
can conclude that the spiral arms are tightly wound. 

Let us find the gravitational potential of the tightly wound 
spiral perturbation (l3^ in the neighborhood of a point 
{Rojfo). Using the properties of WKB approximation, one 
can expand Ei as follows (for more detail see Binney & 
Tremaine 2008) 

Ei(E,4i) (34) 


where 

Ea{Ro) = (35) 

the spiral perturbation (l34l i is reminiscent of a plane wave 
with wavenumber k = keR, where br is the unit vector in 
the radial direction. The potential of a plane wave in a thin 
disk in the framework of the /(E) gravity has been derived 
in Appendix lAl Using equation (IA5b . the potential can be 
expressed as 


$1 


27rGa m§ 

|fc| k"^ +mQ ^ 


(toSI - Uj)k{^a + ha) 

VRa= ^ 

(41) 

2iBk{^a + ha) 

V^ia — ^ 

(42) 

where coefficient functions A and B (the Oort constant of 

rotation) are: 


A = — {mfl — w)^ 

(43) 


(44) 


It is important noting that mathematical form of the Euler 
equation in the weak field limit of /(E) gravity is the same 
as the Newtonian gravity. The only difference between these 
theories appears in the way by which the disk responds to the 
perturbation. On the other hand, the Poisson equation deter¬ 
mines the response or the behavior of the disk to the per¬ 
turbation. Since, for obtaining equations (HTt . the Poisson 
equation has not been used, so we can conclude that they are 
also true in the weak filed limit of /(E) gravity. This is the 
case also for continuity equation (l40t . Therefore, we skip 
the details of the calculations for deriving these equations, 
(l40t and (HU, and refer the reader to chapter 6 of Binney & 
Tremaine 2008. 

Now, we substitute VRa from equation (HU into equation 
(l40t and use equation (l36] l. Also, taking into account that 
ha = uf we find the following dispersion relation 

2 

(mfl — ujY = — 27 rG|fc|Eo q- 

k^ -f rriQ 

We shall derive the local stability criterion from this disper¬ 
sion relation. As expected, in the limit mo —oo and a ^ 1, 
the dispersion relation for WKB modes in Newtonian grav¬ 
ity can easily be recovered. In this paper, we restrict our¬ 
selves to the axisymmetric disturbances, i.e. m = 0. Since 
the right-hand-side of (HU is real, then the disk is stable 
against local axisymmetric disturbances if 0 and un¬ 

stable if < 0. 


(36) 
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It is worthwhile to mention that Toomre’s criteria (l22l i 
and (l 2 ?t have also been derived for axisymmetric perturba¬ 
tions. However, it cannot be concluded strongly that they are 
not applicable for local nonaxisymmetric stability. In fact, 
no general criterion for nonaxisymmetric stability is known 
in Newtonian gravity. 

Now, we consider the general case by adopting non-zero 
speed of sound. In this case, using the dispersion relation 
(Eli, the stability criterion > 0 can be rewritten in the 
dimensionless form 

+ + -‘^\X\+P^ >Q (46) 

Qg 

where 

x = — , /? = ^ (47) 

mo rnaVs 


The schematic behavior of the left-hand side (LHS) of 
equation (|46]) is shown in Figure[T] It is clear from this figure 
that, if the minimum value of the LHS is positive, then the 
condition > 0 will hold for all wavenumbers and conse¬ 
quently throughout the region that we explore, all solutions 
that we will find represent stable modes. In other words, the 
disk is stable against all axisymmetric tightly wound pertur¬ 
bations, if the minimum value of the LHS is positive. It is 
easy to check that the minimum value of the LHS occurs at 
Xmin given by 


Xr, 



3 2 


(48) 


where A is defined as 


A = 


2/3a 



Qg 



2(l-t/32)3Q2 

27/3'-‘a'-‘ 


(49) 
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Fig. 1 Schematic behavior of the LHS of equation ( 146b with re¬ 
spect to X. 


The location of this point varies as the parameters a and /3 
change. Finally the local stability criterion can be expressed 
as 


p ^ 2/3a|Xmm| 

^ X^ . -f (1 -f . -h /32 

min ' I /■‘'-min ' 


(50) 


This is the main result of this section. In fact, equation (fSOl l 
is the generalized version of the Toomre’s local stability cri¬ 
terion for a fluid disk, and should be compared with the stan¬ 
dard case given by (l22l) . Since Xmin contains Qg, criterion 
([soi l is very complicated than the standard Toomres criterion 
(i.e. Qg > 1). However, it is straightforward to show that in 
the limit (3^0 and a ^ 1, the right-hand side (RHS) of 
equation (l50l l equals 1, and the standard Toomre’s criterion 
is recovered. 

It is also instructive to plot the neutral stability curves 
w = 0 for axisymmetric tightly wound waves in a fluid disk. 
These curves specify the boundary of the stable and unsta¬ 
ble modes/waves. To do so, let us introduce a critical wave¬ 
length Xcrit = 2.7:/kcrit- Using this definition and the dis¬ 
persion relation (|45] |. the line separating stable and unstable 
modes is given by 


Qg{y) = a - .2 ( 51 ) 

where y = X/Xcrit and A = 277 jk. We have plotted the 
neutral curve for various values of a and /3 in Figure |2l 
The dashed curves corresponds to Newtonian gravity where 
(a,/3) = (l,0). 

Let us first consider the other curves with /3 = 0, i.e. 
curves corresponding to (q:,/3) = (1 .2,0) and (a,/3) = 
(0.5, 0). In this case one can readily show that the Toomre’s 
criterion (ISOl l is Qg > a. Consequently, if a > 1 (as in 
the case (a, /?) = (1.2,0)) then a lager Qg parameter com¬ 
pared with Newtonian gravity is needed to overcome the lo¬ 
cal gravitational instability. Keeping in mind the physical in¬ 
terpretation of the Jeans instability, this means that the disk 
needs more gas pressure than it would in Newtonian grav¬ 
ity in order to overcome the gravitational collapse. On the 
other hand, if a < 1 then smaller gas pressure is needed for 
preventing the collapse. 

For the general case /3 ^ 0 see the curves (a, /3) = 
(1,0.3) and (a, /3) = (0.8,0.3) in Figure|2] In this case the 
neutral curves corresponding to /? 0 lie below the dashed 

curve. This can not be explained with the above mentioned 
simple interpretation and one should consider it with more 
care. In order to describe this result, for the sake of sim¬ 
plicity and with no loss of generality we assume that a = 1 
and /? 1. We shall see that in reality /3 is small and our 

assumption here is a suitable one. Therefor, the RHS of the 
criterion (fSOl l can be expanded in terms of /3. Consequently, 
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we have 


Qfl > 1 -/?'- + 0(/3") 


(52) 


Thus in the case a = 1, although the gravitational force 
is supposed to be stronger than the Newtonian case, the 
required value of Qg for local stability is always smaller 
than that of the Newtonian gravity. This situation is rem¬ 
iniscent of the difference between Jeans mass in modified 
gravity and Newtonian gravity. In fact in modified theo¬ 
ries of gravity for which the gravitational force is stronger 
than that of Newtonian gravity, the Jea ns mass is smaller 
than the standard case, for ex ample see (IRoshan & Abbassi 
2014 : Capozziello et all 2011 ). 

Therefore, comparing the two cases /3 = 0, a > 1 and 
a = 1, /3 > 0 (note that in both cases the gravitational force 
is stronger than the Newtonian gravitational force), we can 
conclude that in theories where the gravitational force in the 
weak field limit is stronger than that of Newtonian gravity, 
the required value of Toomre’s parameter Qg for preventing 
gravitational collapse can be larger or smaller than the cor¬ 
responding value in Newtonian gravity. In other words, Qg 
in these theories is not necessarily larger than that of New¬ 
tonian gravity. 

In Figure[3]we show the fluid disk dispersion relation for 
growing perturbations for various values of Qg and /3. The 
solid curves correspond to Qg = 0.3 and [3 = 0.1 to 0.5. 
Furthermore, the dashed curves correspond to = 0.5 
and /3 = 0.1 to 0.5. Also, the dotted curves (/3 = 0) are 
the corresponding Newtonian dispersion relations for Qg = 
0.3 and Qg = 0.5. As expected, the growth rates tends to 
increase with decreasing Qg. Also, it is clear form Figure |3 
that the effect of parameter (3 in equation (|45] | is to stretch 
the range of instability to small wavenumber while larger f3 
leads to smaller growth rate. 


3.2 Self-gravitating stellar disk 


Similar to the case of the fluid disk, it is straightforward to 
And the dispersion relation of tightly wound waves in a stel¬ 
lar disk. Since most of the calculations are similar to that of 
the fluid case, we have briefly derived the dispersion relation 
in AppendixlBl The final result is given by equation (IB51 l. In 
the case of axisymmetric perturbations, we may write 


— 27rG|fc|So 


fc2 




OJ 




(53) 


For a given mode of perturbation, if < 0, then the 
mode is unstable. Therefore, one may expect that the system 
is stable against all axisymmetric disturbances if equation 
dSl does not have a solution with < 0. In order to 



Fig. 2 The boundary of stable and unstable axisymmetric WKB 
disturbances in a fluid disk for different values of a and /3. The 
dashed curve corresponds to Newtonian gravity where q = 1 and 

/3 = 0. 



Dimensionless wavenumber -J x 


Fig. 3 Solutions to the dispersion relation ll45t between squared 
growth rate (s'^) and dimensionless wavenumber (kva/n). It has 
been assumed that a = 1. Large p leads to greater instability at 
small wavenumber. 













consider this expectation, we write down equation (l5^ as 
follows 


27rC?|fc|SoaTOQ .to 
{k^ + mQ){K^ — uj^) k ’ 


(54) 


1.0 
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Fig. 4 equation (l55ll in terms of uj^. 


Now, if < 0, i.e. uj = where 7 is a real number, then 
the right-hand side (RHS) of equation (l54l i takes the form 


27rG'|fc|a?7iQEoK ^ 
+ mg) sinh tts ' 



e“^(i+™"^)sinhs'rsinrdr (55) 


where s' = is = 'y/k. Also we have used equation (IB71 i 
for the shape function. The schematic behavior of (l55]) in 
terms of is illustrated in Figure |4] From this figure we 
see that the maximum value occurs at = 0. Therefore, 
if the RHS of equation (l5^ for = 0 is smaller than 1, 
then there is no tightly wound spiral mode with negative 
which satisfies the dispersion relation (l54l i. Consequently, 
the self-gravitating disk will be stable to all WKB perturba¬ 
tions. Therefore, the criterion for stability to all wavelengths 
can be written down as 


27rG|fc|Eoamg ,^ 
{k? + rri^)K^ ’ 


(56) 


this equation can be represented in a more simplified form 
as 


1 (1-e x/o(x)) 

27rGEoa (1 +^ 

where Iq is the modified Bessel function of order 0. Also 
we have used the following identity 



A 

Fig. 5 The solid curve is the exact curve for 5(0, x). and the 
dashed curve is the approximate one, i.e. equation El. 


mx) = -(l-e-^/o(x)) (58) 

X 

Equation (l57l) is the main result of this section. In fact, at 
every location (i?o,(/)o) on the disk, this inequality should 
be satisfied in order to have local stability to all spiral ax- 
isymmetric tightly wound perturbations. As a —1 and 
f3 ^ 0, this criterion approaches to the so-called Toomre’s 
local stability criterion for stellar disk given by (| 2 ^ . again 
as expected. In Newtonian gravity the RHS of (fSTl) is always 
larger than that of metric f{R) gravity. This means that in 
the context of f{R) gravity smaller value of Toomre’s pa¬ 
rameter is required for local stability of rotating disks. As 
we discussed in the previous subsection, this result also is 
the case for a fluid disk. 

Although equation (fSTl i is our final stability criterion for 
every perturbation with wavenumber k, if we And the max¬ 
imum value of the RHS of equation (fSTl) . then we will find 
a criterion which guarantees the stability at every location 
on the disk. Even in Newtonian gravity (where = 0), one 
has to use numerical analysis to find the maximum value of 
the RHS (|52l). Therefore, we need the numeric value of /3 
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to find required criterion. Howe ver, using an exce llent ap¬ 
proximation for reduction factor ( ElmegreenIbOl 1 ). we can 
analytically find the maximum value of RHS of (l57l i. The 
above mentioned approximation is 


(59) 

In the Figure |5l the solid line is the exact value of ^ and the 
dashed line corresponds the approximate value. As it is clear 
from the figure, the deviation between the two functions is 
small. Therefore, with the aid of this approximation, we 
rewrite (|57l) as follows 


CTrK y/x 

27rGEoa {l + P^x){^ + x) 


(60) 


the maximum value of the RHS of this inequality occurs 


-1 - ^2 ^ Vl + 14/32 + ^4 

- w - 

Stability criterion takes the following form 


(61) 


^ ^ QVQay/g{l3) - (1 + / 32 ) 

(5(/3)+5-/f2)(p(^)+5/32-1) ^ ^ 

where p(/3) = \/l + 14/32 + The criterion (l62l i is the 
final result of this section. As expected, the free parameters 
of the theory, i.e. a and mo, appear in the stability criterion. 
If we use the same approximation in Newtonian gravity, then 
we will get the following criterion 


3.14GEo 


> 1 


(63) 


this criterion has a little difference with the exact criterion 
(|23]) . However let us write (l63l l as Qs > 1. As we discussed 
previously, we expect that (3 be small in real situations, so 
the stability criterion can be expanded with respect to (3 as 


Qs> a - a/3^ + + 0{f3^) 


(64) 


this criterion can be compared with the corresponding crite¬ 
rion for a fluid disk (l52] |. The similarity between the stability 
criterion of fluid and stellar disk is obvious. Our discussion 
in the last paragraph of section [TTI considering the parame¬ 
ter Qg for a fluid disk, also holds here. 

In order to make a comparison between the growth rates 
in fluid and stellar disks, we have shown the stellar disk dis¬ 
persion relation (l53T l in Figure |6] The solid curves corre¬ 
spond to Qs = 0.3 and f3 = 0.1 to 0.5. The dashed curves 
correspond to Qs = 0.5 and f3 = 0.1 to 0.5. Also, the 
dotted curves are the corresponding Newtonian dispersion 
relations for Qs = 0.3 and Qs = 0.5 where /3 = 0. As for a 
fluid disk, the growth rates increase with decreasing param¬ 
eter Qs, and the effect of parameter [3 in equation (l5^ is to 


stretch the range of instability to small wavenumber. Also 
larger [3 leads to smaller growth rate. Comparing Figures [3] 
and|6] one can conclude that the the growth rate of perturba¬ 
tion in the stellar disk is smaller than that of the fluid disk. 
In other words, if we consider two patches with the same 
Toomre’s parameter, one in a fluid disk and one in a stellar 
disk, the the growth rate in the patch on the stellar disk will 
be smaller than that of the fluid disk. 


4 discussion and conclusion 


In this paper, the local stability of self-gravitating fluid and 
stellar disks has been investigated in the framework of met¬ 
ric f{R) gravity, for which the function f{R) can be ex¬ 
panded as a power series in terms of Ricci scalar R, see 
equation (|2l). Thus the results of this paper are applicable 
only for these models of f{R). The dispersion relation for 
the propagation of tightly wound spiral waves has been de¬ 
rived analytically for both fluid and stellar disks. The results 
are given by equation (|45]) for fluid disk and (l5^ for the 
stellar disk. 

Also, using the above mentioned dispersion relations, the 
Toomre’s local stability criterion has been derived in the 
framework of metric f{R) gravity. The results are given 
by equations (fSOl) and (l57l i. An important ingredient of the 
present study is its possible application of these criteria in N- 
body simulations of the disk galaxies in the context of f{R) 
gravity. It is worth mentioning that the Toomre’s criterion 
for the fluid disk (fSOl l is complicated than the standard case 
(l22li . On the other hand, the generalized Toomre’s criterion 
for the stellar disk is not so different from its corresponding 
criterion in Newtonian gravity (|2^ . 

As we have discussed, the magnitude of tuq determines 
the magnitude of deviation between local stability criteria 
of f{R) gravity and GR. Therefore, it is necessary to men¬ 
tion that the mass of the effective scalar degree of freedom, 
mg, may depend upon the density of its environment via 
the so -called chameleon mechanism ( Khourv & Weltm^ 
2004 alib i. Obviously, a more careful study is still in order 
to take into account the chameleon mechanism. Studying 
this subject is left as a subject of future study. 
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A Potential of a WKB wave in f{R) gravity 

As we mentioned before, locally, a tightly wound spiral wave can be considered as a plan wave. The reason for this 
similarity is that the pitch angle is very small at every location for tightly wound density waves. It is important to remember 
that the surface density corresponding to a tightly wound spiral wave is given by (l34l) . In order to find the gravitational 
potential of this perturbation, without loss of generality, we choose the x axis to be parallel to k(i?o). Then we use the 
modified Poisson equation (fTSl l to find the gravitational potential. To do so, we guess the solution to (fTSl) as follows 

= A (Al) 


where A and e are arbitrary real constants. The disk is assumed to be thin and so there is no matter outside (z ^ 0) the disk. 
Therefore outside the disk the modified Poisson equation ( fTSl ) can be written down as V"*^<I>i — = 0. Substituting 

the potential (lAll) into this equation, one can easily verify that e = ±k . On the other hand, since matter is located at z = 0, 
derivative of $i with respect to z is not continuous. In order to fix the parameter A, we integrate equation ( |32] | with respect 
to z in the interval z = —^ to z = +^, where ^ is a positive constant, and then let ^ 0. Therefore, the LHS of ( |3^ gives 


f+i 

lim / dz(V^$i 


= lim 




az4 dx'^dz'^ ° dz^ ) 
V dz^ dx^dz ° dz ) 


(A2) 


2|fc|A(P + 


note that we have not written those partial derivatives of $1 which are continuous at z = 0, because their integral with 
respect to z in the above mentioned interval is zero. The same procedure for the RHS of (l32l i gives 


— 47rGamn lim 





[Si(5(z)] 


(A3) 


Taking into account that the tightly wound density wave is Si = and using the properties of the Dirac delta 

function, one can show that the second integral in equation (IA3I) vanishes. Therefore, equation (IA3I) can be simplified as 
follows 




(A4) 


finally, by equating equation (IA3I) to (IA41i . we find the relation between A and as 


27rGa m§ 

|fc| _|_ ^2 “ 


(A5) 


as expected, in the limit mg 00 and a —?► 1, we recover the potential of the tightly wound spiral in Newtonian gravity. 


B Dispersion relation for stellar disk in f{R) gravity 

In order to study the dynamics of a self-gravitating stellar disk in the framework of metric f{R) gravity, we need the 
modified Poisson equation (fTSl l and the collisionless Boltzmann equation (1211) . Using d^TTi . one can readily derive the Euler 
and continuity equations, respectively, as follows 


dvj _ dvj r c)<i) 1 / 2 


(Bl) 






where 


(Jij = ViVj - ViVj 



(B2) 


fvicpx 


(B3) 
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and / is the distribution function. Combining equations (IB lb and (IB2b with (fTSl l and using the same general technique that 
we used in section lTTI to derive the dispersion relation for a fluid disk, we find 


mfl — tu 

VRa = -^ 


(B4) 


where 5^ is the reduction fact or. It is the factor by whic h the response of a stellar disk to an imposed potential is reduced 
below that of a cold disk (see (IBinnev & Tremainel2008h for more detail about this factor and the reason for which it appears 
in the linearized equations), and A has been defined in (l44l i. Also, one can show that equations ( l40b and ( |36] | are also 
applicable for describing a tightly wound density wave in the stellar disk. Therefore, one can find after some algebra that the 
dispersion relation is 


(mVL — loY = — 27rG'|fc| 






(B5) 


where s and x defined as 


uj — mfl 


s = 


X = 


f kan 

V 


(B6) 


As in Newtonian gravity, the difference between the dispersion relation of a fluid disk and a stellar one manifest itself in the 
appearance of the reduction factor in the dispersion relation of the stellar disk. The reduction factor in Newtonian gravity is 
given by (see ( Binnev & Tremain3l2008 ) and references therein) 




1 - 


SinTTS 


3 sin ST sin rdr 


(B7) 


However, a natural question arises here; Does the reduction factor ^ in metric f{R) gravity and Newtonian gravity coin- 
cide? It is worth mentioni ng that in obtaining the reduction factor, the Poisson equation is not used (see Appendix K of 
( Binnev & Tremainel2008 ) for a full derivation of the reduction factor). On the other hand, at least at the level of field equa¬ 
tions, the only difference between Newtonian gravity and the weak field limit of metric f{R) gravity appears in their Poisson 
equation. Therefore, we expect that the reduction factor must be the same in both f{R) and Newtonian gravity. However, 
we need to treat this conclusion more carefully. In fact for obtaining the reduction factor some assumptions have been done, 
and one needs to check their consistency in metric f{R) gravity. This factor has been derived implementing three main 
assumptions. 1) the given perturbation lies in the WKB approximation. 2) The stellar orbit of stars can be described by the 
epicycle approximation in which the orbits are nearly circular. 3 ) The disk is thin and its distribution function is the so-called 
Schwarzschild distribution function (IBinnev & Tremainell2008h . The first assumption is satisfied in this work because we 
are studying tightly wound waves. The main consequence of the second assumpti on which is necessary for deriving the 
reduction factor is the relation between and in the epicycle approximation as ( Binnev & Tremain3l2008 ) 


<^1 I ( dRn{R) \ 

cr^ 2 \ d\nR J 


(B8) 


this equation is derived only using the kinematic properties of the nearly circular motion and does not depend on the grav¬ 
itational theory. Thus this relation can also be used in modified gravity theories including f{R) gravity. Finally, for the 
third assumption it is just needed to note that, as we showed in subsection 12.21 the mathematical form of the collisionless 
Boltzmann equation in f{R) and Newtonian gravity is the same. Consequently, one can readily verify that the Schwarzschild 
distribution function is also a solution for (l2Tl i. Therefore, we can be sure the reduction factor (IB7b is also true in metric 
/(i?) gravity. 
















